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ABSTRACT 



Context. Previous analyses of lithium abundances in main sequence and red giant stars have revealed the action of mixing mecha- 
nisms other than convection in stellar interiors. Beryllium abundances in stars with Li abundance determinations can offer valuable 
complementary information on the nature of these mechanisms. 

Aims. Our aim is to derive Be abundances along the whole evolutionary sequence of an open cluster We focus on the well-studied 
open cluster IC 4651. These Be abundances are used with previously determined Li abundances, in the same sample stars, to investi- 
gate the mixing mechanisms in a range of stellar masses and evolutionary stages. 

Methods. Atmospheric parameters were adopted from a previous abundance analysis by the same authors. New Be abundances have 
been determined from high-resolution, high signal-to-noise UVES spectra using spectrum synthesis and model atmospheres. The 
careful synthetic modeling of the Be lines region is used to calculate reliable abundances in rapidly rotating stars. The observed 
behavior of Be and Li is compared to theoretical predictions from stellar models including rotation-induced mixing, internal gravity 
waves, atomic diffusion, and thermohaline mixing. 

Results. Beryllium is detected in all the main sequence and turn-off sample stars, both slow- and fast-rotating stars, including the 
Li-dip stars, but is not detected in the red giants. Confirming previous results, we find that the Li dip is also a Be dip, although the 
depletion of Be is more modest than for Li in the corresponding effective temperature range. For post-main-sequence stars, the Be 
dilution starts earlier within the Hertzsprung gap than expected from classical predictions, as does the Li dilution. A clear dispersion 
in the Be abundances is also observed. Theoretical stellar models including the hydrodynamical transport processes mentioned above 
are able to reproduce all the observed features well. These results show a good theoretical understanding of the Li and Be behavior 
along the color-magnitude diagram of this intermediate-age cluster for stars more massive than 1.2 Mq. 

Key words, stars: abundances - stars: evolution - stars: rotation - open clusters and associations: individual: IC 4651 



1. Introduction 

The classical theory of stellar evolution, which only allows for 
mixing in stellar convective layers, fails to explain nearly all Li 
abundance patterns observed so far at the surface of low-mass 
stars. Indeed, in contradiction to the classical predictions, the 
vast majority of F- and early G-type stars (including the Sun) 
deplete, often quite severely, their surface Li abundance during 
their main-sequence lifetime, revealing transport processes be- 
yond convection. A-type stars are also concerned, although in 
these objects the signatures of non-standard Li depletion, occur- 
ring in their interior during the main sequence, appear at their 
surface only when they cross the Hertzsprung gap and become 

subgiants. 

Lithium depletion is seen bot h in field (e.g ., iHerbigiil965|; 
lDuncanll98l1:lPasquini et al.ll994tlObre et al.ll 999; Chen et al] 
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Jones et al.|[T997l iPasquiniet al1l200U ISestito & Randichll2005 
and references therein). As cluster stars have well-defined 
masses and, most important, share the same age and initial chem- 
ical composition, they are ideal targets for investigating mixing 
processes as a function of stellar mass and evolutionary status. 
They have therefore received considerable attention in the liter- 
ature. Observations of the Li evolution along the color magni- 
tude diagram (CMD) of open clusters ha ve revealed a number of 
interesting features, srich as the Li-dip (IWallerstein et al.lll965l ; 
iBoesgaard & Tripiccolll986[). the lack of abun dance trends with 
age among A-type stars liurkhart & Couprvl [T998. 2000), the 
presence of a spread o f Li among solar-type stars in relativel y 
old clusters like M67 dPasquini et al.lll997l: Hones et al.lll999l) . 
and the Li variat ions with mass and e volutionary status as in the 
case of IC 465 1 (iPasquini et al.ll2004 . 

Many different physical mechanisms have been pro- 
posed to exp l ain these observat ions; atomi c diffusion 
(iMichaudI [19861: iMichaud et al.ll2004l) . mass loss dHobbs et al.l 
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Table 1. Observational data of the sample stars. 
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Fig. 1. Color magnitude diagram of IC 4651. Sample stars are 
shown as filled circles. The ubxy photometry is from iMeibomI 
(I2OOOI) and obtained by means of the WEBDA database. 



[19891 ^wenson &Faulkneii Il992l). rotation-induced mix- 
g JCharbon nel et al J 1 1 9921: Delivan nis & Pinsonneaulj 



19971 iTalon & Charbonnell 119981: IPalacios et al] I200I 



Theado & Vau clait< l2003h7 mi xing by internal gravity waves 



dGarcia Lope" & SpruiT fl99ll: iMontalban & SchatzmanI l2000t 
lYoung et al.llToOSi) . or a combination of some of these processes. 
So far, only the hydrodynamic stellar models that combine 
the effects of meridional circulation, shear turbulence, internal 
gravity waves, and atomic diffusion have been able to account 
for Li observations over a broad range of stellar masse s, ages, 
and evolutionary status jTalon & Charbonnell I2003L 120051: 
ICharbonnel & Talonll2005L l2008h . These models, which use a 
prescrip tion for the wave excitation that reproduces the solar-p 
modes jGoldreich et al.l 1 19941) . explain other observational 
constraints concomitantly, like the behavior of C, N, O, and Be 
across the Li dip in young cluster s, such as the Hyades, or th e 
internal rotation profile of the Sun (ICharbonnel & Taionll2005h . 

Additional clues can help ascertain the nature and descrip- 
tion of transport mechanisms of chemicals and angular momen- 
tum inside low-mass stars. In particular the simultaneous deter- 
mination of Li and Be abundances in the same stars strength- 
ens the constraints by performing a stellar tomography. Such a 
powerful diagnosis is possible because Li and Be burn at differ- 
ent temperatures (~ 2.5x10'' K for Li and ~ 3.5x10'' K for Be), 
whi ch correspond to differ ent depths in the stellar interior (see 
e.g. lDelivannis et al.ll2000l) . 

Although the simultaneous study of Li and Be in the 
same stars is a powerful tool. Be has been determined in 
far fewer stars than Li has. In particular, despite the advan- 
tages of cluster stars mentioned above . Be has been in vesti- 
gated i n only a few open cluster s (e.g. iBoesgaard et al.l 



l2004al: iGarcia Lopez et al.l 1 1 9951: iRandich et alj 12002 " 



1977 



2007 



and references therein) , while field star s have been more ex- 
tensively studied (e.g. IBoesgaardI [197 6: Stephens et alj 119971 : 
ISantos et alj|2004l: IBoesgaard & Krugler Holleki,2009l and ref- 
erences therein). 

Boesgaard and collaborators discovered a Be dip and a cor- 
relati on between Li and Be depletion in mid F-stars in the 
field dPelivannis et al. 19981) and also in the Hyades, Coma, 



Star 


V 


(b - y) 


obs. date 


'-exp (^) 


S/N 


E3 


12.05 


0.38 


10.may.2001 


2 X 2100 


90 


E5 


12.82 


0.33 


09.may.2001 


2 X 3000 


90 


E7 


14.17 


0.40 


26.mar.2000 


4 X 3600 


60 


E14 


13.12 


0.34 


09.may.2001 


2 X 3300 


95 


E15 


13.52 


0.38 


10.apr.2001 


2 X 4500 


85 


E19 


12.24 


0.42 


08.apr.2001 


2 X 2100 


85 


E34 


13.42 


0.36 


ll.may.2001 


2 X 4500 


90 


E45 


14.18 


0.44 


29.mar.2000 


3 X 3600 


40 


E56 


12.18 


0.38 


08.jun.2001 


2x2100 


85 


E64 


13.72 


0.37 


06.jun.2001 


2 X 4500 


100 


E79 


13.55 


0.37 


09.apr.2001 


2 X 4500 


90 


E86 


13.74 


0.39 


06.jun.2001 


2 X 4500 


100 


E95 


11.95 


0.52 


08.jun.2001 


2 X 2100 


80 




LZ.ia 


U. j4 


08.jun.2001 


2 X 2100 


on 


T1228 


12.94 


0.35 


08.jun.2001 


2 X 3300 


90 


E25 


12.65 


0.37 


07.jun.2001 


2 X 2100 


85 


T2105 


13.95 


0.43 


Randich et al. (2002) 




E8 


10.67 


0.66 


10.may.2001 


1500 


50 


E12 


10.32 


0.64 


15.apr.2001 


2 X 1200 


75 


E60 


10.87 


0.67 


09.apr.2001 


2 X 1800 


70 


E98 


10.88 


0.66 


ll.may.2001 


2 X 1800 


65 


T812 


11.06 


0.67 


07.jun.2001 


2 X 1800 


70 



The Stromgren photometry is from Meibom (2000) . The S/N, per 
resolution element (4 pixels), was measured in the Be region. The 
five stars grouped at the end of this and the following tables are the 
red giants. 



l2004al) . They also show that for cooler G stars there is little 
or no Be depletion, even i n Li-depleted objects up to ~ 5500 
K jBoesgaard et al.l l2004bl an d references therein). For even 
cooler stars, below ~ 5500 K, ISantos et all (|2004 show a Be 
decline correlated with decreasin g temperature. These findings 
were confirmed bv lRandich et al.| (|2002, 2007i) . who investigated 
Li and Be in late F- and early G-type stars in some open clusters 
over a relatively broad age range (IC 2391 - 50 Myr, NGC 2516 
- 150 Myr, Hyades - 600 Myr, IC 4651 - 1.7 Gyr, and M67 - 
4.5 Gyr). 

To improve our understanding of the mixing processes, 
adding new Be data covering a large fraction of the color mag- 
nitude diagram of well-studied open clusters is important. The 
ambitious aim is to compare the observations with models for 
stars over as wide as possible a range of mass and evolutionary 
status. In this context, the open cluster IC 4651 is a very good 
test c ase. It has well-determined membership dMeibom et al.l 
I2OO2I) . and precisely kno wn metallicity and chemical com- 
position ([Fe/H1 = -1-0 .11: iPasquini et"al] 12004 ICarretta etaP 
12004 iPaceetal I I2OO8I) . Age determinations for this cluster vary 



and Praesepe clusters jBoesgaard & King|2002l : lB"oesgaard et alj 



between 1.2+ 2 Gy r dBiazzo et alj |2007|) and 1.7+0.15 Gyr 
dMeibom et alj|2002l) . A color magnitude diagram of IC 4651 
is shown in Fig.[T] 

Important is that the Li evolution along the CMP of this 
cluster has already been investigated in detail by iPasquini et al.l 
d2004 . Well-determined patterns have^been found along its 
evolutionary sequence (see Fig. 3 by Pasquini et al. 2004-). In 
a sequence of increasing stellar mass, the cluster indeed has 
an Li "plateau" around solar mass stars, followed by a well- 
determined Li dip, after which the maximum Li content is ob- 
served, at meteoritic value. A sudden Li drop follows along the 
Hertzsprung gap, possibly indicating the onset of dilution (the 
so-called first dredge-up). Some of the giants show noticeable 
differences in Li abundance while sharing other stellar char- 
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acteristics. It has been sho wn t hat models of rotatin g stars by 
ICharbonnel & TalonI (Il999l) and lPalacios et al.l (l2003h were able 
to reproduce the Li behavior both of main sequence stars lying 
on the blue side of the dip and of evolved stars. 

In this work we present Be a bundances for the same stars an- 
alyzed by iPasquini et alj (l2004h . This is the first time Be abun- 
dances have been determined for stars along the whole evolu- 
tionary sequence of an open cluster, including solar-type, Li-dip, 
turn-off, subgiant, and red giant stars. With these new results, the 
mixing processes in different stellar masses can be investigated 
in greater detail. 



2. Observations and stellar parameters 

Our sample is co mposed of 22 stars. Tw enty one of them have 
been analyzed by iPasquini et al.l (|2004|) . To these we a dd on e 
star of this same cluster analyzed by iRandich et al.l (l2002h . 
star T2105. The sample was selected to cover a large inter- 
val of the CMP from the Stromgren photometry obtained by 
MeibomI (l2000h. The interested reader is referred to the work of 



Pasquini et al. ( 120041) for more details on the sample selection. 



Spectra were obtained using UVES, the U ltraviolet and 
Visual Echelle Spectrograph jPekker et alj|2006h fed by UT2 of 
the VLT. UVES is a cross-dispersed echelle spectrograph able 
to obtain spectra from the atmospheric cut-off at 300 nm to ~ 
1100 nm. The reduction was conducted with the ESO UVES 
pipeline within MIDAS. The spectra have typical signal-to-noise 
(S/N) between 40 and 100 and R ~ 45 000. The log book of 
the observations is given in Table [T] Stars with names s tarting 
with an "E" follow the numbering system of'Eggenll'lQTl!). Stars 
with names starting wi t h a "T " follow the numbering system of 
lAnthonv-Twarog et all (Il988l) . 

The atmos pheric parameters for the sample stars determined 
by IPasquini et al. (2004.) were adopted. A first estimate of the ef- 
fective temperature (Teg) w as calculated using the calibrations of 
lAlonso et alj ( ll996lll999[) . A first estimate of the surface gravity 
was obtained using this Teff, a mass of 1.80 Mq for the cluster 
turn-o ff, and the distance modulus determined bv lMeibom et al.l 
(I2OO2I) . (m-M)o = 10.03. Starting with these values, both the Teff 
and the microturbulence velocity (^) were further constrained 
using Fe i lines, eliminating trends with the excitation potential 
and the equivalent widths, respectively. The ionization equilib- 
rium, on the other hand, was not adopted to constrain the gravity. 
While the other parameters were being varied, gravity was also 
varied by at most 0.30 dex to obtain as good an equilibrium as 
possible without departing much from the log g given by the star 
position at the CMD. 

Using this method. jPasquini et al.l(l2004h had some difficulty 
constraining the parameters of the fast-rotating hot stars around 
the turn-off. The above-mentioned approach resulted in higher 
metallicities for these stars than the ones determined for the 
slow-rotating stars. This effect was likely caused by the blending 
of the Fe i fines with other neighboring lines due to the rotational 
broadening. Thus, the measured equivalent widths were greater 
than they should be. As in these cases the spectroscopic param- 
eters are likely not reliable, atmospheric parameters for the fast- 
rotating stars were calculated adopting the photometric temper- 
atures. The spectroscopic temperatures were adopted for the re- 
maining stars. The atmospheric parameters of the stars are given 
in Table [2] The mean m etallicity of the cluster as determined 
by IPasquini et al.l ( |2004|) was adopted for all the stars, [Fe/H] - 
H-0.11. 




3129.5 3130 



3130.5 
X(k) 



3131 



3131.5 



Fig. 2. Fit to the region of the Be lines in the main sequence 
star E7. The solid line represents the observed spectrum, the 
dotted line a synthetic one with solar abundance ratios, and the 
dashed line a synthetic spectrum with [0/Fe] = H-0.25 and [C/Fe] 
= +0.10. 



3. Be abundances 

3.1. Synthetic spectra 

For late-type stars. Be abundances can only be determined using 
the Ben -S - ^Pq resonance lines at 3131.065 A and 3130.420 
A. This near-UV region is extremely crowded with atomic and 
molecular lines, some of them still lacking proper identification. 
Thus, determination of Be abundances needs to be conducted 
with spectrum synthesis taking all the blending nearby features 
into account. 

To derive the Be abundances, sy nthetic spectra were cal- 
culated with the codes described in ICoelho et al.l (l2005b and 
the grids of model atmosphere s without overshooting ca lculated 
with the ATLAS9 program dCastelU & Kuruczl l2003b . These 
models assume local thermodynamic equilibrium, plane-parallel 

geometry, and hydrostatic equilibrium. 

Th e line list is the same as used in ISmiljanic et al.l (l2008l 
2009bl). In short, the molecular line list is described in 



Coelho et al. ("2005^ and the atomic line list is the one com- 
piledby lPrimas et all (1 19971) . Log gf of -0. 168 and -0.468 were 
adopted for the Be lines at 3131.065 A and 3130.420 A, respec- 
tively. The line list includes an Fe i line in 3131.043 A, with log 
gf - -2.517 and x - 2.85 eV, affecting the blue wing of the 
Be 3131 line. The parameters of this line were constrained us- 
ing several stars o f different parameters and metallicities (see 
iPrimas et al.ll 19971 for details). The proper identification of this 
line, however, remains controversial in the literature, and other 
choices have been made by differe nt authors. As discussed in 
more detail bv lRandich et al.l (|2007|) . this choice has no signifi- 
cant effect on the calculation of Be abundances of stars with Tgff 
> 5400 K. 



3.2. The Sun and the slow-rotating stars 

That some of the sample stars are fast rotators adds further dif- 
ficulty to the analysis of the Be spectral region. Thus, we con- 
ducted a careful analysis, first modeling the slow-rotating stars 
and conducting a few tests to make certain any possible system- 
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Table 2. Stellar parameters and abundances of Li and Be 



Star 




Spec. 


log g 




1' sini 


A(Be) 


A(Li) 




(K) 


or Phot. 


(dex) 


(km s-') 


(km s-i) 


(dex) 


(dex) 


E3 


6550 


Phot. 


3.90 


2.10 


29.9 


1.21 


1.64 


E5 


6930 


Phot. 


4.20 


2.00 


23.9 


1.36 


3.47 


E7 


6300 


Spec. 


4.30 


1.10 


2.1 


1.11 


2.83 


E14 


6860 


Phot. 


4.20 


1.90 


34.1 


0.76 


<2.07 


E15 


6540 


Phot. 


4.20 


1.70 


10.0 


0.64 


<1.93 


E19 


6280 


Phot. 


3.90 


2.10 


32.1 


0.41 


<1.66 


E34 


6640 


Phot. 


4.20 


1.90 


24.2 


0.11 


<1.92 


E45 


6350 


Spec. 


4.30 


1.10 


4.2 


1.08 


2.80 


E56 


6520 


Phot. 


3.90 


2.10 


27.8 


1.21 


<2.16 


E64 


6650 


Spec. 


4.30 


1.70 


9.2 


0.76 


2.24 


E79 


6620 


Spec. 


4.30 


1.70 


21.8 


0.51 


<1.91 


E86 


6600 


Spec. 


4.30 


1.70 


14.0 


0.96 


2.20 


E95 


5800 


Spec. 


3.50 


1.70 


12.0 


0.61 


2.19 


E99 


6830 


Phot. 


4.00 


2.00 


28.1 


1.21 


<2.38 


T1228 


6770 


Phot. 


4.20 


1.70 


5.7 


1.21 


3.18 


E25 


6680 


Phot. 


4.00 


2.00 


21.3 


1.11 


3.29 


T2105 


6110 


Phot. 


4.44 


1.10 


8.0 


1.11 


2.82 


E8 


4900 


Spec. 


2.70 


1.30 


0.1 




<0.09 


E12 


5000 


Spec. 


2.70 


1.50 


1.1 


<-0.70 


<0.35 


E60 


4900 


Spec. 


2.90 


1.40 


1.9 


<-0.70 


0.42 


E98 


4900 


Spec. 


3.00 


1.40 


1.0 


<-0.50 


0.72 


T812 


5000 


Spec. 


3.00 


1.60 


0.5 


<-0.70 


<0.38 



atic effect affecting the modeling of the fast-rotating stars would 
be under control. 

As the first step, with the line list described above, we fitted 
the Be lines in the solar UVESQ spectrum, adopting the param- 
eters: Teff = 5777 K, log g = 4.44, and ^ = 1.00 km s^'. An 
abundance of A(Be) = 1.10 was obtained. T his abundance is in 
excelle nt agreement with the one found by IChmielewski et all 
(Il975h . A(Be) = 1.15, usually adopted as the reference photo- 
spheric solar abundance. 

The low photospheric solar Li abundance, depleted by a 
factor of ~ 160 with respect to the meteoritic valu e, has been 
known for a long time (see e.g. iMiiller et al.lll975l) . The pho- 
tospheric Be abundance and its difference with respect to the 
meteoritic one have instead been a matter of debate. All anal- 
yses of the solar photospheric abundance obtain log(Be/H) ~ 
1.15-1.25 (IChmielewski e"tanfT97l iGarcia Lopez et alj|l995t 



iRandich et all 20021) , a factor o f 2 less than w hat is measured 
in meteorites, log(Be/H) = 1.41 (|Lodders"2003'). This issue now 
seems settled to a large extent; the measured photospheric value 
is lower than the real one due to our inability of properly ac- 
counting all continuum opacity sources in the region around 313 
nm, where the two Be ii res onance lines are located, when calcu- 
lating model atmospheres (iBalachandran & Bellll 19981 lAsplundl 

The second step was to fit the two main sequence slow ro- 
tator stars E7 and £45. The cluster wa s found to have a so- 
lar abundance pattern bv iPasquini et al.l (|2004). Thus, the solar 
abundance ratio^ were adopted when calculating a first set of 
synthetic spectra. However, as shown in Fig. [2l when adopting 



' The spectrum is available for download at the ESO website: 
www.eso.org/observing/dfo/quality/UVES/pipeline/solar_spectrum 
.html 

^ The abundances recommended bv lGrevesse & Sauva l( fl998l) were 
adopted for all elements except for oxygen, for which we adopted the 
abundance suitable f or the ID atmospheric models recommended by 
lAllende Prieto et alj j200lh . A(0) = 8.77, and for Be for which the 
abundance derived above from the solar spectrum was adopted. 



a solar abundance pattern, an important molecular feature near 
the 3130 A Be line is not properly fitted. A better fit is obtained 
when [O/Fe] is increased by +0.25 dex and [C/Fe] by -i-0. 10 dex, 
in both E7 and E45. These apparent overabundances are likely 
not to be real but caused by shortcomings in the analysis, such 
as the neglect of effects caused by departures from the local ther- 
modynamic equilibrium and/or by problems with the molecular 
and atomic data in this region. 

Nevertheless, it is important to notice that properly fitting 
the neighboring molecular features in the slow-rotating stars is 
an important step toward analyzing the fast-rotating ones. In the 
last case, all these lines are blended into a single feature. As all 
the stars in the cluster are expected to have the same chemical 
composition (except for the giants), one may expect that the en- 
hanced abundances of C and O are also needed to properly fit 
the spectra of the fast-rotating stars. We stress that this increase 
in the abundances of C and O has no effect on the calculated Be 
abundances of the slow-rotating stars. Thus, when fitting the fast 
rotators, C and O abundances were kept with these increased val- 
ues, while the other elements were kept with a solar ratio. Only 
the Be abundance was varied from star to star The derived Be 
abundances are listed in Table |2] 

3.3. The fast-rotating stars 

Before fitting the true fast-rotating stars, we conducted a test 
to better understand the effects of the broadening in the de- 
rived abundances. We artificially broadened the spectra of the 
two slow-rotating stars, E7 and £45, by convolving them with 
Gaussian functions. In Fig. [3] we show the spectrum of the star 
£7 broadened with four different Gaussian functions. The Be 
abundance was then calculated for these broadened spectra. 

In the four cases shown in Fig. |3] an abundance of A(Be) 
- 1.16 dex was derived for £7, while A(Be) = 1.11 dex was 
derived in the original unbroadened spectrum. Applying similar 
broadening to £45, we obtain A(Be) =1.15 for all the broad- 
ened cases, while A(Be) = 1.08 dex was obtained for the origi- 
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Fig. 3. Fit to the region of the Be lines in the main sequence star E7. The four panels show a sequence of larger broadening applied 
to the observed spectrum. In all panels the best-fitting synthetic spectrum has A(Be) = 1.16 dex. 



nal observed spectrum. Thus, we conclude that the abundances 
in a broadened spectrum are slightly overestimated. However, 
this small amount is well within the uncertainties (see below). 
Therefore, we are confident that it is possible to derive mean- 
ingful abundances for the fast-rotating stars with our method of 
analysis. 

Beryllium abundances were then derived for all the fast- 
rotating main-sequence and turn-off stars. The stars were found 
to have a broad range in abundances, from 1.36 dex (E5) to 
0.11 dex (E34). In Fig. H] we plot the spectra of these two stars, 
which have similar v sini, to show that the difference in Be abun- 
dance obtained in the analysis is real. In the lower panel of this 
same figure we show the result of the division of these spec- 
tra. The structures left in the division have central wavelengths 
of 3130.41 and 3131.03 A, arguing that the difference is really 
caused by the different Be abundances. Similarly, we compare 
in Fig.|5]the observed spectrum of star E34 and the spectrum of 
E7 artificially broadened to a v sini similar to that of E34. It is 
clear from this figure that star E34 has a smaller Be abundance 
than E7. This comparison again shows that the difference in Be 
abundances detected in the stars is real. 

Finally, we point out that Be was detected in all the sam- 
ple of main sequence and turn-off stars. Most of the stars show 
signs of Be depletion, i.e.. Be lower than the meteoritic value, 
but still have a detectable amount of Be in their photospheres. 
In Fig. |6]we show the observed spectra of E19 and E34, the 
two stars with smaller Be content. Also shown are a best-fit syn- 
thetic spectrum together with a spectrum calculated with no Be 
at all. As clearly seen, the spectra lacking Be show a marked 
difference when compared to the observed ones. This compari- 
son argues that our Be abundances are actual detections and not 
upper-limits. The derived Be abundances are listed in Table |2] 



E34 vs E5 




3128 3130 3132 3134 

A (A) 



Fig. 4. (a) Comparison between the observed spectra of stars E5 
(dashed line) and E34 (solid line). The stars have similar rota- 
tional velocity and different Be abundances, (b) The result of the 
division of the two spectra in the upper panel. The two resid- 
ual features are located at 3130.41 and 3131.03 A. These wave- 
lengths are compatible with the difference being caused by dif- 
ferent Be abundances. 

3.4. The red giants 

Fitting synthetic spectra in the Be region for the red giants 
proved to be even more challenging than fitting the fast-rotating 
stars. From this exercise it became clear that the line list used 
to fit the main sequence and turn-off stars is not fully adequate 
for red giants. Problems with the wavelengths and depths of the 
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Fig. 5. Comparison between the observed spectrum of star E34 
(solid line) and the spectrum of star E7 (dashed line) artificially 
broadened to a v sini similar to the one of E34. It is clear the 
lower Be abundance of star £34 and that a higher abundance 
could be easily detected in our analysis. 



Fig. 7. Fitting the Be lines in the red giant E98. The crosses 
represent the observed spectrum, the solid line the best fit, the 
dashed lines represent variations of + 0.20 dex from the best fit, 
and the dotted lines represent a spectrum without Be and one 
withA(Be) =1.21. 
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Fig. 6. (a) Comparison between the observed spectrum of star 
E34 (solid line), the best synthetic fit (dotted line), and a 
synthetic spectrum calculated without Be (dashed line), (b) 
Comparison between the observed spectrum of star El 9, the best 
synthetic fit, and a synthetic spectrum calculated without Be. 
These two comparisons argue for both stars having detectable 
amounts of Be. 



lines become clear. More work on the transitions affecting the 
near UV region is required. 

In particular, the Fe i line adopted as the blend affecting the 
3131 A Be line appears to be too strong here and with a wrong 
wavelength when compared to the observed spectra (see Fig.|7]i. 
For the dwarf stars, as discussed before, this blend has no signif- 
ica nt effect. For metal-poor stars, such as the subgiants analyzed 
by iGarcia Perez & PrimasI ( l2006l) , the blending feature is also 
not significant. For our metal-rich cool giants, however, it pre- 
vents reliable abundances from being derived from this Be line. 
In this case, we relied solely on the 3130.4 A feature. 



The 3130.4 A line, however, is very weak (indicating that 
Be is strongly depleted or even absent in giants) and strongly 
blended. In Fig. [7] an example of fitting the Be region in a red 
giant is given. The Be line only affects the wing of a very strong 
and heavily blended feature. The best fit is found when the wing 
of this feature in the synthetic spectrum is brought to the level 
of the wing in the observed spectrum. The figure suggests that 
an abundance cannot be derived in this case with an uncertainty 
better than +0.20 dex. 

One should also note that the abundances in this case are 
strongly dependent on broadening, blending features (their in- 
clusion and also on the right abundance for the given element), 
and on the determination of the continuum. Thus, one should be 
very careful when interpreting the abundances derived for the 
red giants. In this case we consider that only an upper limit on 
the Be abundance can be determined. In addition, the difficul- 
ties discussed above mean that one should be very careful not to 
push the comparison between dwarfs and giants too far The Be 
upper limits for the red giants are listed in Table |2l 

Of all these issues, the level of the continuum is of particu- 
lar concern. Different solutions could bring the wing of the line 
up or down, changing the best fit. Finding the correct level of 
the continuum in the near UV region is at best a very difficult 
task. However, by adopting very similar solutions for all the stars 
that have very similar atmospheric parameters, we believe that at 
least a reliable differential scale can be obtained. 

Regarding the abundances of other elements, it is expected 
that the first dredge-up would change C and N abundances of 
the stars on the red giant branch (RGB). This near UV region is 
affected by many CN, NH, and CH lines, so this change in abun- 
dances could be an important issue. We calculated a synthetic 
spectrum for star E60 by adopting typical post dredge-up values 
for these elements, [C/ Fe] = -0.15 and [N/Fe] = +0.40 (see e.g. 
ISmilianic et al. I l2009al) . A comparison with the best fit obtained 
using unaltered abundances showed only very weak effects, well 
within the +0.20 uncertainty of the fit. 

The Be abundances of all sample stars are shown as a func- 
tion of the V magnitude in Fig. [8] For comparison, we also show 
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Fig. 8. The Be abundances of the sample stars, calculated in this 
work, as a function of the V magnitude. Abundance determina- 
tions are shown as solid circles and upper limits as open trian- 
gles. The Be dip is clearly seen. 



Fi g. 9. The Li abundan ces of the sample stars, calculated 
in iPasquini et al.l ( |2004|) . as a function of the V magnitude. 
Abundance determinations are shown as solid circles and upper 
limits as open triangles. 



Table 3. The uncertainties of the Be abundance introduced by 
the uncertainties of the atmospheric parameters and by uncer- 
tainty of the fitting itself related to the level of the continuum 
and the S/N. 



Star 


O" Tcff 


O-logg 




O" [Fe/H] 




O" total 


E7 


±0.03 


±0.10 


±0.00 


±0.00 


±0.05 


±0.12 


E86 


±0.05 


±0.10 


±0.05 


±0.05 


±0.10 


±0.17 


E60 


±0.15 


±0.05 


±0.03 


±0.05 


±0.20 


±0.26 



a plot of the Li abundances determined bv lPasquini et al.l ( |2004|) 
as a function of V in Fig.|9] The evolutionary stages of the stars 
are indicated in the figures. It is clearly seen that the stars of the 
Li dip also define a clear Be dip. 

3.5. Uncertainties and comparison witli tine literature 

The determination of Be abundances is affected by the uncer- 
tainties of the atmospheric parameters and by the uncertainty in 
the determination of the pseudo-continuum during the spectrum 
fitting. To conduct the calculations of the errors, we adopted the 
same uncertainties of the atmospheric parameters determined by 
IPasquini et all (|2004 . crxeff = + 100 K, crjogg = + 0.30 dex, cr^ = 
+ 0.30 km s"', and cr^Pe/H] = ± 0.10 dex. 

To estimate the effect of the atmospheric parameters, we can 
change each one by its own error, keeping the other ones with the 
original adopted values, and recalculate the abundances. Thus, 
we measure what effect the variation of one parameter has on 
the abundances. We assume the uncertainties introduced by each 
parameter to be independent of the other ones. The uncertainty 
due to the continuum was determined by estimating the sensitiv- 
ity of the derived Be abundance on the adopted continuum level. 
This uncertainty is mostly related to the S/N of the spectrum. 
These effects are listed in Table [3] The calculations were done 
for three representative stars, E7 as an example of slow-rotating 
star, £86 as an example of a fast-rotating star, and E60 as an 
example of a giant star. 



The Li and Be abundances derived for the star E5, which 
seem to be slightly above the solar values (see Table |2|, seem to 
suggest that the uncertainty on its temperature is slightly higher 
than for the other stars. For these hot stars, the effective temper- 
ature is also important for determining Be abundances. We thus 
believe that its Li and Be abundances are slightly overestimated 
and should be interpret ed with care. 

As mentioned in ^ iRandich et al.l (I2OO2I) has analyzed three 
main sequence stars of IC 465 1 and determined their Be abun- 
dances. Two of their stars, E7 and E45, a re also analyzed here . 



The atmospheric parameters derived by iRandich et al. 
are different from the ones adopted here. iRandich et al 



(2002) 



(2002) 



calculated temperatures of 6061 and 6016 K, while we adopted 
6300 and 6350 K for stars E7 and E45, respectivel y. A gravity 
of log g - 4.44 was adopted by lRandich et al.l (l2002h for all stars 
they analyzed, whil e a value of log g = 4 .30 was calculated for 
both E7 and E45 in IPasquini et al.l (|2004|) and adopted here. In 
spite of th ese differences, the Be abundances of the two stars as 
derived by IRandich et al.l ( 20021) agree w ith ours within the un- 
Randich et al ] j2002l) derived A(Be) = 1.11 

1.11 



certainties; that is 
and 1.16 for stars E7 and E45 
and 1.08, respectively. 



while we obtained A(Be) 



4. Comparison with model predictions for rotating 
stars 

4.1. Input physics of the hydrodynamical stellar models 
4.L1 . Hot side of the Li dip 



In IPasquini et al.l (l2004l) the Li abundance patterns in IC 465 1 
were compared with model predictions for rotating stars by 
Charbonnel & TalonI (Il999l) and Palacios et al. (2003[). These 
models include atomic diffusion, as well as the transport 
of angular momentum and chemicals by meridional ci rcula- 
tion a n d shear turbul e nce f ollowi ng the formalism by ZahnI 



(119921) . iTalon & ZahnI (1 19971) . and lMaeder & ZahnI (1 19981) (for 
more details and refe rences see the original papers as well as 
iDecressin et al.l 120091) . They are adequate for stars that lie on 
or originate in the hot side of the dip and inside which angu- 
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lar momentum transport by intern al gravity waves emitted b y 
convective envelopes is inefficient ( iTalon & Charbonnell2003l) . 
These models were found to explain the corresponding Li data 
for both mai n sequence and gian t stars in the intermediate-age 
IC 465 1 (see iPasquini et al.ll20()4h . as they do for younger clus- 
ters like the Hyades, Com a Berenices, and Praesepe (see e.g. 
ICharbonnel & Talonll2008l) . 

Meanwhile new stellar models within the same mass range 
(i.e., for stars more massive than 1.2 Mq) were computed for 
other purposes with the stellar evolution code STAREVOL V3. 1 
(Charbonnel & Lagarde, in preparation) using the same pre- 
scriptions for the transport of angular momentum and chemi- 
cals by meridional circulation and shear turbulence as in the 
above-mentioned papers. Atomic diffusion is also included in 
the form of gravitational settling, as well as that related to ther - 
mal gradients, using the formulation by iPaquette et al.l (Il986l) . 
As far as the i nput physics is concerned, t he di fferences with 
the m odels by ICharbonnel & TalonI (Il999h and iPalacios et al.l 
(l2003b are the fol lowing: (1) the reference solar composition is 
the one derived bylAsplund et al.' ("2006"), except for Ne for which 
the recommendations by Cunha et al. ( 2006?) are followed; (2) 
new OPAL opacity tables are used according to the chemical 
mixture; (3) thermohaline mixing is accounted for following 
ICharbonnel & ZahnI (|2007|) . a process that only becomes effi- 
cient when the stars reach the so-called luminosity bump on the 
RGB (see SJi. 

Results of these recent computations for initial stellar masses 
of 1.3, 1.35, 1.4, 1.5, 1.8, and 2 Mo for Z = 0.014 are presented. 
In view of the remaining uncertainty on the age of IC 4651, 
we provide in Table Irrelevant model characteristics at 1.2 and 
1 .7 Gyr for the stars that are still on the main sequence at that 
time. In all cases various initial rotation velocities on the zero- 
age main sequence were considered. (The 1.3, 1.35, 1.4, and 
1.5 Mq models are computed for Vzams - 50, 80, and 110 km 
s"'; the 1.8 M© models are computed for Vzams = 80, 110, and 
180 km s-\ and the 2 Mq models for Wzams = 1 10, 180, and 250 
km s The efficiency of magnetic braking applied to stars ar- 
ri ving on t he ma in sequence is determined according to th e work 
of lGaigel ( Il993h foll owing the bri&ing l aw by Kawaler ( Il988h 
(for more details see lPalacios et al.|[2003h . In practice, magnetic 
braking only has to be applied to the 1.3, 1.35, and 1.4 M© stars 
ffiat have thicker convective envelopes than more massive stars. 
As in previous works, the aim was to reproduce the rotation ve- 
locity distribution observed at the age of the Hyades. Theoretical 
rotation velocities at the age of IC 4651 are given in Table|4|and 
span the velocity range of our sample stars. In all cases the sur- 
face rotation velocity also changes because of secular evolution, 
and it is reduced to only a few km s ' when the models are on 
the subgiant branch due to the radius expansion. 

4.1.2. Cool side of the Li dip 

For less massive main-sequence stars lying on the cool side of 
the Li dip that have deeper convective envelopes, internal gravity 
waves are efficiently generated and participate in the extraction 
of internal angular momentum, along with meridional circula- 
tion and shear turbulence (Talon & Charbonnel 2 003). In this 
case, t he ap propriate models are those by iTalon & Charbonrieil 
(I2005h and Icharbonnel & TalonI (|2005|) . which take all these 
effects into account. Unf ortunately the 1.0 M© models of 
ICharbonnel & TalonI ( l2005h have an effective temperature at the 
age of the IC 4651 (i.e., 5597 K), lower than our cooler main 
sequence stars, and they cannot be used for comparison with 
the present sample. Obviously abundance determinations for 



such relatively cool objects would be extremely valuable for 
further tests of the models. Only the 1.2 M© model computed 
byljalon & Charbonnel (2005) is available for comparison pur- 
poses (see 94. 2b . It has an initial rotation velocity on the zero 
age main sequence of 50 km s"', and surface spin-down follows 
Kawaler (1988) as explained in § |4XT] 

4.1 .3. Initial values for Li and Be 

In the following we compare the predictions of the hydrodynam- 
ical models described above with the observations of both Li and 
Be in IC 465 1 . Results are presented separately for the main se- 
quence and turnoff stars ( 94.21 ) and for more evolved stars ( 94.3t . 
In the following figures, we use the V magnitudes to divide the 
stars according to their evolutionary stages; stars fainter than V 
- 12.6 are at the main sequence, stars between 12.6 and 1 1.8 are 
at the upper turn-off or at the subgiant branch, and stars brighter 
than 11.8 are giants. 

For the initial values of A(Li) and A(Be) used in our 
computations we ch ose the met eoritic abundan ces, i.e., 3.25 
(lAsplund et all |2006|) and 1.41 (iLoddersI l2003h . respectively. 
This Li value is actually lower than the highest Li abundance 
determined among our sample stars, which induces an artifi- 
cial shift within the comparison between predictions and data. 
For this reason we increase, within the figures, the theoreti- 
cal Li points by 0.16 dex. This corresponds to assuming that 
the highest Li abundance derived among our sample stars (E5, 
A(Li)=3.47 dex) is the initial abundance of the sample stars. 

4.2. Comparisons with observed features on the main 
sequence 

Figures [TOj and [TT] compare the predictions of the stellar mod- 
els described in 94.1! to the Li and Be observed abundances of 
the main-sequence and turn-off stars of our IC 465 1 sample as 
a function of the effective temperature. The models provide an 
estimate of the expected dispersion in surface abundances due 
to a dispersion in initial rotation velocities. On the cool side of 
the dip, only one theoretical point is available in the temperature 
range consid ered in the present study (1 .2 Mq, Vzams = 50 km 
s ' model of lTalon & Charbonnell2005b . 

It is possible to note on these graphs (see also Fig. [T2b 
that the observational behavior of the Be abundances follows 
that of the Li abundances closely. Th e Li dip is also a Be dip , 
confirming previous lite rature results dBoesgaard & Kin3l2002l ; 
iBoesgaard et al.ll2004ah . In a narrow range in effective temper- 
ature centered on ~ 6650 K, stars present both an Li and a Be 
depletion. Important is that Li, which is more fragile (it burns 
indeed through proton-captures at ~ 2.5x10^ K compared to ~ 
3.5x10^ K for Be), is more strongly depleted within the dip. In 
most of the stars lying in this region, only upper limits could be 
derived for Li, while Be is always detected. 

The predictions of the models described in § 14.1! success- 
fully account for the Li patterns as a function of effective tem- 
perature, namely, the presence of a large fraction of stars with Li 
upper limits between ~ 6500 and 6800 K (i.e., within the dip) 
and the higher Li values on both the cool and hot sides of the 
dip (see Fig.lTOl). They also reproduce the Be behavior remark- 
ably well in the whole effective temperature range considered 
(see Fig.nri). The success of the models is confirmed when plot- 
ting Be versus Li as shown in Fig.[T2l As noted in § 4.1.3, the 
theoretical Li values have been shifted by 0.16 dex to account 
for the difference between the meteoritic Li and the highest Li 
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Table 4. Rotation velocity, effective temperature, and surface abundances of Li and Be predicted for the stellar models that are still 
on the main sequence at the age of IC 4651 (models by Charbonnel & Lagarde, in preparation). See text for more details. 



M 


VzAMS 


V/ 








A(Li) 


A(Be) 


(Mo) 


(km .?"') 


(km 




(K) 
















1.2 Gyr 


1.7 Gyr 


1.2 Gyr 


1.7 Gyr 


1.2 Gyr 


1.7 Gyr 


1.2 Gyr 


1.7 Gyr 


1.3 


50 


19.4 


16.9 


6532 


6507 


2.05 


1.51 


1.02 


0.77 




80 


28.3 


24.4 


6535 


6528 


1.24 


0.37 


0.62 


0.22 




110 


36.2 


30.9 


6520 


6483 


0.47 


-0.59 


0.2 


-0.29 


1.35 


50 


30.3 


27.6 


6687 


6663 


2.69 


2.33 


1.26 


1.10 




80 


28.6 


24.5 


6683 


6670 


1.16 


0.24 


0.57 


0.15 




110 


36.6 


31.5 


6685 


6676 


0.43 


-0.70 


0.20 


-0.35 


1.4 


50 


42.7 


41.4 


6834 


6762 


3.24 


3.22 


1.41 


1.41 




80 


67.8 


67.8 


6824 


6764 


3.20 


3.10 


1.41 


1.39 




110 


96 


95.2 


6826 


6781 


3.12 


2.94 


1.39 


1.33 


1.5 


50 


40.5 


37.5 


7094 


6857 


3.23 


3.18 


1.41 


1.40 




80 


66.5 


62.7 


7174 


6948 


3.13 


2.94 


1.40 


1.34 




110 


92.2 


87.8 


7136 


6900 


2.97 


2.66 


1.35 


1.21 




6600 
Teff 




Fig. 10. Li abundances in IC 465 1 main sequence and turnoff 
stars fainter than V ~12.6 (black points and triangles for ac- 
tual determinations and upper limits, respectively). Open circles, 
squares, and triangles show model predictions for initial rotation 
velocities of 50, 80, and 110 km s"', respectively. For compari- 
son purposes the theoretical Li points are increased by 0.16 dex 
to account for the difference between the meteoritic value and 
the highest Li abundance determined in our sample stars. See 
text for more details. 



abundance determined in our sample stars; however, what really 
matters for constraining the models is the slope of the trend be- 
tween Li and Be. As can be seen in Fig. [12] the theoretical and 
observational slopes are found to agree very nicely. Linear fits 
of the observational values, accounting for errors in both axes, 
have slopes of -1-0.58 + 0.09, if all detections and upper lim- 
its are considered, of -1-0.45 ± 0.05 if we only exclude star £34 
(the one with the smaller Be abundance), and of -1-0.35 + 0.08 
if we consider only the detections. For the theoretical values 
we calculate a slope of H-0.45 considering all points. If we re- 
strict the theoretical values to the ones with A(Li) > 2.0, which 
is the range where the observational detections are, a slope of 
-1-0.33 is obtained. These values agree very well with the ob- 
served ones. Similar correlated depletions of Li and Be have 



Fig. 11. Same as FiglTO|for observed and theoretical Be abun- 
dances. 



also been observed in field and cluste r stars by DeUvannis et al.l 
(ll998h . lBoesgaard et al.1 (l2004bh . and lRandich et al.l (l2007h . 

On the cool side of the dip (i.e., below ~ 6400 K), we 
find no Be dispersion thus confirming previotis resii lts by 
iBoesgaard et"!!] (l2004bl) and iRandich et all (l2002l |2007|) . This 
provides further constraints to models of rotating low-mass (i.e., 
below ~ 1 .25 Mq), stars including the transport of angular mo- 
mentum by gravity waves that should assume a range in initial 
rotation velocities (Charbonnel et al. in preparationfl However, 
we can already note that models of solar-mass stars computed 
with a large range of initial angular momentum (V„„ from 15 to 
110 km s ') and including the efifects of internal gravity waves 



' iRandich et al.l j2002h compared the Li and Be abundances de- 
termined in early G-type stars in old clusters, among which 
IC 4651 with predictions of models including gravity waves by 
iMontalban & SchatzmanI l l2000h . In these models, rotation-induced pro- 
cesses are ignored, and waves are invoked as the only source of mixing 
for chemicals in the stellar radiative zone. This is very different from 
the view proposed by Charbonnel & Talon. In this case indeed gravity 
waves transport angular momentum, together with meridional circula- 
tion and turbulence, but the transport of chemicals is essentially caused 
by rotation-induced mixing. See the original papers for more detailed 
information. 
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12 3 4 

A(Li) 

Fig. 12. Li and Be abundances in IC 465 1 main sequence and 
turnofF stars fainter than V ~12.6 (black points and triangles for 
Li actual determinations and upper limits, respectively). Open 
circles, squares, and triangles show model predictions for initial 
rotation velocities of 50, 80, and 110 km s"', respectively. For 
comparison purposes the theoretical Li points are increased by 
0.16 dex to account for the difference between the meteoritic 
value and the highest Li abundance determined in our sample 
stars. See text for more details. 




8000 6000 1000 8000 6000 4000 
I'ell 'I'cfl 

Fig. 13. Li and Be abundances in IC 465 1 subgiant stars with 
Teff between ~ 6830 and 5800 K and V between 12.6 and 11.8. 
Theoretical predictions for the surface Li and Be evolution are 
shown for the 1.8 M© star (see text for more details). The solid 
line is for the classical case, while the other lines present the 
results for the rotating models computed with different initial 
rotation velocities (80, 110, and 180 km s"', respectively, for 
dashed, dotted, and dot-dashed lines). The Li tracks are shifted 
by 0.16 dex as explained in the text. 



predict very small Li dis persion at the age of IC 4651 (see 
ICharbonnel & Talonll2005h . and thus even more modest Be dis- 
persion. This is caused by internal gravity waves that dominate 
the transport of angular momentum in stars lying on the red side 
of the dip and enforce quasi solid-body rotation within the stel- 
lar interior As a result, the surface Li and Be depletions are ex- 
pected to almost be independent of the initial angular momentum 
distribution (i.e., of the initial rotation velocity), implying very 
low surface abundance variations from star to star On the other 
hand, in more massive (i.e., higher than ~ 1.25 Mq) stars the 
efficiency of internal gravity waves dramatically drops and in- 
ternal differential rotation is expected to be maintained along the 
stellar life under the effects of meridional circulation and turbu- 
lence. Consistently, variations in the initial angular momentum 
from star to star lead to more Li and Be dispersion, as required 
by the observations. 

4.3. Comparisons with features observed in subgiant and 
giant stars 

Let us now focus on the post-main-sequence evolution. 
Figures [13] and [14] compare the Li and Be abundances of our 
subgiant and giant sample stars, as a function of the effective 
temperature (which now depicts evolutionary stage), with pre- 
dictions of our most massive models computed by assuming dif- 
ferent initial rotation velocities. Standard (i.e., non-rotating) pre- 
dictions are also shown. The subgiant stars (with Tetf between ~ 
6830 and 5800 K) are compared with the 1.8 model pre- 
dictions, while the relevant tracks for the giants (with Tea at ~ 
4900-5000 K) are the 2Mo ones. 

In the classical (non-rotating) case, the surface depletion 
of both Li and Be begins relatively late in Ten within the 
Hertzsprung gap (i.e., at ~ 5500 K). This depicts dilution of the 
external layers with nuclear-processed material when the con- 



vective stellar envelope deepens in mass. However and as shown 
in Fig. [13] actual dilution starts earlier than predicted classi- 
cally, as already revealed by L i data both in ope n clusters like 
IC 46 51 jPasquini et al.1 1200 4') and NGC 3680 ('Pasquini et ajJ 
20011). but a l so in field stars (jA ls chuler 1975; Bro wn et al. 1981 



Lebre et al.l 119991: iDo Nascimento et al l 120001: iPalacios et al.l 
120031) . As can be seen in Fig. \T3\ the Be depletion also occurs 
earlier than predicted in the classical case. This simply means 
that the Li- and Be-free regions are larger in real main sequence 
stars than classically predicted. 

This favors the complete hydrodynamical models, since 
rotation-induced mixing acting on the main sequence enlarges 
the Li- and Be-free regions inside the stars, implying earlier (i.e., 
at higher Teu inside the Hertzsprung gap) Li and Be depletion at 
the surface of subgiant stars, as well a s lower Li and Be va lues 
after dredge-up in giant stars (see also iPalacios et al.ll2003L for 
model comparison with Li data for field stars). The behavior of 
both elements predicted in the rotating case nicely matches the 
data for the subgiant and giant stars in IC 4651, and a disper- 
sion in initial rotation velocity explains the observed abundance 
dispersion along the evolutionary sequence well. 

Lithium abund ances of the red gian ts of IC 4651 are dis- 
cussed in detail bv iPasquini et al.l (|2004|) . As can be seen from 
Table [2] Li was detected in only two giants, although all five of 
them have very similar atmospheric parameters. To explain this, 
Pasquini et al. (2004) suggests that the Li observed in these stars 
might have been produced in their interior after the first dredge- 
up. We can see, however, that this is not necessary, since in the 
corresponding stellar mass domain, differences in initial rotation 
velocities are expected to lead naturally to a significant disper- 
sion (~ 1.6 dex for the range in initial rotation assumed here) 
in post-dredge-up Li values, accounting nicely for the observed 
data. 
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Fig. 14. Li and Be abundances in IC 465 1 evolved stars brighter 
than V ~ 12.6. Predictions are shown for the 2 Mq model. The 
solid line is for the classical case, while the other lines present 
the results for the rotating models computed with different initial 
rotation velocities (110, 180, and 250 km s"' for dotted, and dot- 
dashed, and dot-long dashed lines, respectively). The Li tracks 
are shifted by 0. 16 dex as explained in the text. 

iPasquini et al.l (|2004|) also discusses the issue of the evolu- 
tionary status of the giant stars in IC 465 1 . Their proper classi- 
fication is based only on their position at the CMD and is not 
straightforward. These stars could indeed be either on the first 
ascent RGB or at the clump. As mentioned in ^4.1.11 the new 
models presented here are computed to inc lude thermohaline 
mixing following ICharbonnel & Zahnl(l2007l) . This process only 
becomes efficient when the stars reach the bump on the RGB. It 
then leads to an additional decrease in both the Li and Be surface 
abundance at Teft ~ 4250 K, as can be seen in Fig. [14] Clump 
stars are thus expected to have lower Li and Be surface abun- 
dances than their RGB counterparts. However, for a 2 Mq star 
the effect is relatively modest (see also Charbonnel & Lagarde, 
in preparation), and the theoretical predictions cannot help us 
disentangUng the evolutionary status of the giants. 

5. Conclusions 

Beryllium abundances were calculated for twenty-one main se- 
quence, turn-off, subgiant, and giant stars belonging to the open 
cluster IC 465 1 . This is the first time that such an analysis has 
been carried out along the whole evolutionary sequence of an 
open cluster. 

The Be abundances are found to clo sely follow the behav- 
ior of the Li abundances as determined bv lPasquini et al.l (|2004|) . 
The coolest main sequence stars observed present no Be abun- 
dance dispersion, and a well-defined Be dip overlaps the well- 
documented Li dip for F-type stars. The Be abundances of post- 
main-sequence stars confirm that first-dredge up dilution starts 
earlier than the expected classically and present a significant dis- 
persion. This is consistent with the Li behavior in evolved stars 
in IC 4651, in NGC 3680, and in the field. 

Our observations are compared with theoretical predictions. 
For the main sequence and turnoff stars lying on the hot side 
of the Li-dip, as well as for the evolved stars, we compared the 
observations with results of new hydrodynamical stellar mod- 
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els calculated with STAREVOL V3.1 (Charbonnel & Lagarde, 
in preparation). These models include the effects of atomic dif- 
fusion, meridional circulation, shear turbulence, and thermoha- 
line mixing. (Wave-induced transport is negligible in the case 
of these stars that have very narrow convective envelope while 
on the main sequence.) For the main sequence stars lying on 
the cool side of the Li dip, we compa red the observations with 
the 1.2 Mq model of ICharbonnel & "faloa (.2005i) which also in- 
cludes the transport of angular momentum by internal gravity 
waves. 

The models reproduce all the Li and Be features very nicely 
along the CMD of IC 465 1 . The dispersion of Li and Be abun- 
dances on the blue side of the dip and in evolved stars is very 
well explained when accounting for a dispersion in the initial 
values of the stellar rotational velocity as observed in young 
open clusters. On the other hand, the lack of abundance disper- 
sion in lower mass stars is expected to be caused by the impact 
of internal gravity waves. 

The success in explaining the Li and Be abundances along 
the whole evolutionary sequence of IC 465 1 is very encourag- 
ing. It shows that important steps have been taken towards the 
proper understanding of the physical mechanisms acting during 
the stellar evolution. 
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